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Abstract 

The algebraic-geometric approach is extended to study solutions of iV-component systems 
associated with the energy dependent Schrodinger operators having potentials with poles in the 
spectral parameter, in connection with Hamiltonian flows on nonlinear subvariaties of Jacobi 
varieties. The systems under study include the shallow water equation and Dym type equation. 
The classes of solutions are described in terms of theta-functions and their singular limits by 
using new parameterizations. A qualitative description of real valued solutions is provided. 
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1 Introduction. 



The quasi-periodic solutions of most classical integrable PDEs can be obtained using the inverse 
scattering transform (1ST) (see, for example, Ablowitz and Segur [1981], Newell [1985] and Ablowitz 
and Clarkson [1991]). This is done by establishing a connection with an isospectral eigenvalue 
problem for an associated Schrodinger operator. 

The solution of nonlinear evolution equations using algebraic-geometric techniques was initially 
developed to handle A^-phase wave trains. This approach can be summarized as follows. By using 
the trace formula, families of quasi-periodic and soliton solutions are associated with Hamiltonian 
flows on finite dimensional phase spaces. These flows are described by using so called /i-variable 
representations leading to an Abel-Jacobi mapping which include holomorphic and, in some cases, 
meromorphic differentials (see amongst others, Ablowitz and Ma [1981], Dubrovin [1981], Ercolani 
[1989] and Alber and Alber [1985]). Then the mapping is inverted in terms of Riemann theta- 
functions and their singular limits. Many well-known nonlinear equations such as KdV, sine- 
Gordon, focusing and defocusing nonlinear Schrodinger equations, which describe a wide variety 
of important phenomena in physics, optics, biology and engineering, were studied by using this 
approach. 

Recently special attention was given to the shallow water (SW) equation derived in Cammasa 
and Holm [1993] in the context of the Hamiltonian structure, 

U t + 3UU X = U xxt + 2U X U XX + UU XXX - 2kU x , (1.1) 

and the Dym type equation (see Cewen [1990], Hunter and Zheng [1994] and Alber et al. [1994, 
1995]) 

U xx t + 2U X U XX + UU XXX — 2kU x = 0, k = const. (1.2) 
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Camassa and Holm [1993] described classes of n-peakon soliton-type solutions for an integrable 
(SW) equation ( |l.l|) . In particular, they obtained a system of completely integrable Hamiltonian 
equations for the locations of the "peaks" of the solution, the points at which its spatial derivative 
changes sign. In other words, each peakon solution can be associated with a mechanical system of 
moving particles. Calogero [1995] and Calogero and Francoise [1996] further extended the class of 
mechanical systems of this type. 

The problem of describing complex traveling wave and quasi-periodic solutions of the equations 



(1.1) and ( |1.2j ) can be reduced to solving finite-dimensional Hamiltonian systems on symmetric 
products of hyperelliptic curves. Namely, according to Alber et al [1994,1995,1999], such solutions 
can be represented in the following form 

U(x, t) = m H + fi g -M, 

where g is a positive integer, M is a constant and the evolution of the variables \i is given by the 
equations 

9 ..fcj... f k = l,...,g-2, 
i=l 2 V R (Vi) ( dx k = g. 

Here i?(/i) is a polynomial of degree 2g + 2 (for shallow water equation (|1.1|)) or 2g + 1 (for the 
Dym type equation (1.2)). Also M = for the Dym type equation. 



In contrast to the finite-dimensional reductions of such equations as KdV and sine-Gordon 
equations, system Ql.3| ) contains a meromorphic differential, also the number of holomorphic differ- 
entials is less than the genus g of the corresponding hyperelliptic curve: W 2 = R((i). This implies 
that the problem of inversion ( |1.3| ) can not be solved in terms of meromorphic functions of x and t. 
Examples of such equations arise in several problems of mechanics. These were considered in Van- 
haecke [1995] and Abenda and Fedorov [1999], where a connection was established with the flows 
on nonlinear subvarieties of hyperelliptic Jacobian varieties, so-called strata. In Alber et al. [1997] 
a whole class of A-component systems with poles was shown to be integrable by reducing them 
to similar nonstandard inversion problems which contained meromorphic differentials. Therefore 
A-component systems can be over deter mined, implying that the genus of the spectral curve can 
be higher then the number of /^-variables. 

iV-component systems can be briefly described as follows. For the KdV equation, the spectral 
parameter appears linearly in the potential of the corresponding Schrodinger equation: V = u — A, 
in the context of the 1ST method. In contrast, Antonowicz and Fordy [1987a, b, 1988, 1989] and 
Antonowicz, Fordy and Liu [1991] investigated potentials with poles in the spectral parameter 
for what they refer to as energy dependent Schrodinger operators connected to certain systems 
of evolution equations. Specifically, they obtained multi-Hamiltonian structures for A-component 
integrable systems of equations related to the following isospectral eigenvalue problem: 

/ d 2 V\ 

L Ha? + *r =0 ' (L4) 

M N 
3=0 3=0 

where the kj are constants and the Vj(x, t) are functions of the variable x, the parameter t and the 
spectral parameter A is complex. This includes the coupled KdV and Dym systems. 
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In Alber et al [1994, 1995, 2000a], the presence of a pole in the potential was shown to be 
necessary for the existence of weak billiard solutions of nonlinear equations. Billiard solutions of 
nonlinear PDE's have been related to finite-dimensional integrable dynamical systems with reflec- 
tions including ellipsoidal Birkhoff billiards. It turned out that the existence of billiard solutions 
and the presence of monodromy effects is a specific feature of the whole class of iV-component 
systems with poles (see Alber et al. [1997]). 

Quasi-periodic solutions of the Dym equation were studied in Dmitrieva [1993a] and Novikov 
[1999] by using a connection with KdV equation and introducing additional phase functions. Soliton 
solutions of the Dym type equation were studied in Dmitrieva [1993b]. Periodic solutions of the 
shallow water equation were discussed in McKean and Constantin [1999]. Beals et al. [1999, 
2000] used Stieltjes theorem on continued fractions and the classical momemt problem for studing 
multi-peakon solutions of the (SW) equation. 

The main goal of this paper is to describe explicit formulae in terms of theta-functions and their 
singular limits for the solutions to the shallow water equation (1.1). We also explain the role of the 
mysterious phase functions used by Dmitrieva [1993a] when studying SW equation and equation 
of the Dym type. This phase is present for the whole class of WKI hierarchy, (see Wadati et al. 
[1979]). 

In the present paper the traveling wave, soliton, peakon, cuspon and quasi-periodic solutions are 
considered. Usually in the case of integrable evolution equations quasi-periodic flows are liberalized 
on the Jacobi varieties. In this paper we show that in the case of N-component systems with 
poles the x- and i-flows take place on nonlinear subvarieties (strata) of generalized (noncompact) 
Jacobians. For this reason, the term "liberalization" is no longer applicable here. This makes the 
above nonlinear equations quite different from such well known equations as KdV, sine-Gordon 
and nonlinear Schrodinger equations. For the sake of clarity, in this paper we start with solutions 
related to (hyper)elliptic curves of at most genus 2. The case of arbitrary genus is only notationally 
more complicated and we provide complete formulae. In addition, we give a complete classification 
of real bounded solutions U(x, t) in the above cases and provide corresponding plots. 

Notice that the complex geometry of the traveling wave solutions, cusp and peakon solutions was 
previously studied in Alber et al [1994, 1995, 2000] in connection with geodesic flows with reflections 
on Riemannian manifolds and in Li and Olver [1998] from the point of view of singularity analysis. 



The Contents of the Paper. In Section 2 we demonstrate the main difference between the 
nonlinear SW equation and the KdV equation from the point of view of the algebraic-geometric 
approach by obtaining traveling wave solutions as a result of inverting elliptic integrals of the second 
and third kind. Here we express the amplitude U and the phase X as meromorphic functions of 
an auxiliary variable parameterizing the elliptic curve. 

In Section 3 we apply different singular limits to the problem of inversion resulting in formulae 
for different periodic and solitary solutions. In particular, peakon solutions are obtained as limits 
of the traveling wave solutions and are related to various singularizations of the elliptic curves. 

Section 4 provides explicit expressions in terms of theta-functions for the stationary quasi- 
periodic solutions. This is done by using new complex parametrizations on the related associated 
hyperelliptic curve of genus 2. 

In Section 5 we find time-dependent solutions by integrating and inverting equations ( |1.3[ ) 
in the genus 2 case. We show that these equations can be extended to a standard Abel-Jacobi 
mapping of a symmetric product of the hyperelliptic curve to its generalized Jacobian. The original 
system (|1.3j ) then defines a mapping onto a 2-dimensional nonlinear stratum of the Jacobian, a 
generalized theta-divisor, where the dynamics actually takes place. By fixing t in the expression 
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for the solution in terms of theta- functions, we then recover the stationary quasi-periodic solutions 
obtained in Section 4. 

Section 6 contains qualitative analysis of real bounded solutions for the case when the Weier- 
strass points of the spectral curve are real. 

In a forthcoming paper we will consider different singular limits of the quasi-periodic solutions 
when the spectral curve becomes singular and its arithmetic genus drops to zero. The solutions 
are then expressed in terms of purely exponential tau-functions and, in the real bounded case, they 
describe, in particular, a quasi-periodic train of peakons tending to a periodic one at infinity. 



2 Traveling Wave Solutions. 

The main difference between the nonlinear SW or Dym equations and the KdV equation from 
the point of view of the algebraic-geometric approach can be demonstrated already on the level of 
traveling wave solutions. 

The traveling wave solutions of the KdV equation are obtained by inversion of an elliptic integral 
of the first kind, i.e., the integral of a holomorphic differential, which results in a meromorphic 
function. 

In contrast to this, after substituting U(x, t) = X(x — ct) into the Dym type equation ( |1.2D and 
using a simple transformation (see Alber et al [1995]), we arrive at the problem of inversion of the 
following Abelian integral of the second kind 




dX = f — \ ° d\ = x-ct = X, (2.1) 



k(X - ai)(X - a 2 ) J Uo ^^R 3 (X) 

defined on the elliptic curve £ = {W 2 = Rs(X)}, where 

R 3 (X) = (X - ai )(X - a 2 )(X - c) 

and Uo being a constant. Here we suppose that the roots of Rs(X) are distinct. The differential 
in ( |2.1[ ) has a double pole at infinity A = oo and a double zero at the Weierstrass point A = c 
on £. It follows that the complex inverse function U(X) must have two independent periods on 
C, the periods of the differential along two nontrivial homology cycles on £. On the other hand, 
U(X) blows up only at X = oo. There are no meromorphic functions with such properties (see e.g. 
Markushevich [1977]). Moreover, because of the double zero of the differential, the solution U(X) 
has moving branch points of the form 

U -c = 0({X -X ) 2/3 ). (2.2) 

One can show that U(X) has infinitely many such branch points. 

To deal with this we describe the complex function U(X) in a new parametric form by intro- 
ducing a new complex variable u which gives a parameterization of the curve £ : 

dX 

u. (2.3) 



/a 2/R^A) 
Then, according to the theory of elliptic functions, 



d 2 

X(u) = p(u) + A = — — r log An (z) + const, z = 2iriu/uj3 (2.4) 
du z 
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where p(u) is the elliptic Weierstrass function with periods 2uj\ and 2u>3 depending on the coeffi- 
cients of i?3(A) and A = (oi + ai + c)/3. The #11(21) is the quasi-periodic Riemann theta- function 
which vanishes at the points of the period lattice 
A = {z = 2iriZ + BZ}, B = (see, e.g. Dubrovin [1981]): 

e u {z) = exp Qb(M + 1/2) 2 + (M + l/2)(z + mU . 
Here we use two variables u, z because of different normalizations of (quasi)-periods of p and 9. 



The curve £ can be identified with the factor C/A. Now the integral ( |2.1| ) can be transformed as 
follows 



/ (p(u) + A — c) du = + const, no = const, (2-5) 



which yields 



JkX + const = C(u) + (A - c)u = ^- log 9 n (z) + A'z, A' = const, (2.6) 
where £(u) is the Riemann zeta-function with parameters 771 and 773 such that for any u £ C 
C(u)= p(u)du, C(« + 2wi) = C(w) + 2»7i, C(w + 2^3) = C(n) + 2r/ 3 . 



The constants 2r/i, 2773 are interpreted as the periods of the differential 

Xd\/2y / R 3 (\) = X(u) du. 

Then, in view of (2.6), 

VkX(u + 2wi) = y/JzX(u) + 2r}\ + 2(A - c)wi, 
\/kX(u + 2w 3 ) = v^X(u) + 2773 + 2(A - 0)^3. 



Thus we have expressed the amplitude U = X and the phase X in terms of the auxiliary complex 
variable u. 

In what follows we study real solutions U(X) = X(u(X)) which correspond to the case of all roots 
of the polynomial i?3(A) being real. We choose the half-period oj\ to be real and uj^ to be purely 



imaginary. The existence of branch points of U(X) (see (2.2)) implies that the real solutions may 



have cusps. To demonstrate this we consider two different cases. 

Case 1: k > 0, a± < 02 < c < 00. According to the theory of elliptic functions, parameterization 
(HJD yields that 

A(cji) = c, \(u> 3 ) = at, A(ui + lj s ) = a,2, A(0) = 00. (2.7) 

Along the real axis $tu and the line {u = 103 + u'\u' € R}, A (it) is real valued. Since we are 
interested in nonsingular traveling wave solutions and X(u) has a pole at the origin, let us consider 
X(u)3 + u') as a function of u'. It is periodic with period 2^3 and takes values in the interval [a\, 0,2]- 



The differentials in (2A) and ( |2.5| ) do not have poles or zero's on the line {u = UJ3 + u'}. Hence, 
X(u>s+u') and the inverse u'(X) are monotonic functions. As a result, due to (|2.6|), the composition 
function U(X) = A(ll>3 + u'(X)) is a regular periodic function with period (2773 + 2(A — c)ui 3 ) / y/H, 
and its graph is merely a distortion of that of A(u>3 + v!) . 
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Case 2: k > 0, a\ < c < a 2 < oo. Here we have 

= a 2 , X(u; 3 ) = a 1 , X(ui + a> 3 ) = c, A(0) = cxd. (2.8) 

Again, A(u) is real valued along the real axis 3iit and the line {u = + ii'|u' G R}. Since X(u) 
has a pole along the real axis, we again consider only the real function X(loi + u'), u' G R which 
is 2cj3-periodic and now takes values in the interval [c, 02]. In this case, the differential Q2.5|) has a 
double zero at if = u>i + (u' = <^i), where the derivative dX/du has a simple zero. As a result, 
the derivative dU/dX = dX/du ■ du/dX blows up for the corresponding value of X, which implies 
that the graph of the function U(X) has cusps with periodicity (2% + 2(A — c)u>3)/t/k . This 
phenomenon was first detected in Alber et al [1994,1995]. 



SW Equation. Traveling wave solutions for the shallow water equation have a similar description. 
Substituting U(x,t) = fi(x — ct) into ( |1.1| ), we obtain an Abelian integral of the third kind 



U J . J C . . '-dfi= [ U /i / ° dfi = x-ct = X, (2.9) 

IU V «(M - a l)(^ - a 2)(V ~ a 3) Ju y/Ky/Rdjl) 

defined on the even order elliptic curve 

£ = {W 2 = Ri(fJt)}, Ri(^) = (M - ai)(/i - a 2 )(/i - a 3 )(/i - c), 



where the roots of R4 are supposed to be distinct. The differential in ( |2.9D has a pair of simple 
poles at the infinite points oo_, oo + on E. Then the complex inverse function U(X) must have three 
independent periods on C: the periods of the differential along two nontrivial homology cycles on 
£ and along a homology zero cycle around one of the infinite points. A meromorphic function with 
such a property does not exist. 

Let A, B be canonically conjugated cycles on £ and u) be the normalized holomorphic differential 

da 

uj - 



U^RjJxj' 



where the multiplier II is chosen from the condition £. uj = 2ni. Introduce a new variable z 
parameterizing £ as follows 

'''' z, (2.10) 



where 01,02 and c denote the corresponding Weierstrass (branch) points on £. Then we get the 
expression 

d , e[d](z-q/2) 

* ) ->-di> a ' m, + w (2 - u> 



(id), 



dfi 



T A Jc n^R 4 {») 

Thus fj,(z) is an elliptic function with periods 2iri, B = <f B u. The integral ( |2.9| ) takes the form 



v^X = /%(*) - c) dz = (p - c)z - log ° ni * q/ ^ + const. (2.12) 
Jz e n {z + q/2) 
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As a result, we have expressed the function U = fi and its argument X in terms of z. 

Since the differential in (|2.9| ) has a double zero at the Weierstrass point [i = c, similar to the 
case of the Dym equation, the inversion of ( |2.9| ) yields U{X) with infinitely many branch points of 
the form fl2.2p . This results in real solutions having cusps. 

Now we describe 2 different types of real traveling wave solutions of assuming that k > 
and all roots of R/±{[i) are real. In this case the period B is real as well. 



Case 1. a\ < a2 < < c. According to the parameterization (|2.10|) , we have 



q e 



n 



a ■> 



Ol 



Mo) 



, MB/2) 



B 



M(<z/2) = Z7(- 9 /2) 



a 3 n^A^/x) 
/i(7ri + B/2) = a 2 , /i(vri) = a 3 . 



oo, 



The function ^(z) is real along the real axis 5iu and the line £ = {z = ni + z'\z' E E}, but it is 
finite only in the second case. The function fj,(-Ki + z') changes periodically between a 2 and 03 with 
period B. Along the line i the differential in ( |2.S| ) has no zeros, therefore X{iti + z') and the inverse 
real function z'{X) are strictly monotonic functions. According to ( 2.12| ), the composition function 
U(X) = fj,(ni + z'(X)) is a regular periodic function with real period 



(p - c)B - log 



MB-g/2) 
MB + <?/2)J ' 



Such a function is shown in figure la. 



Case 2. a\ < a 2 < c < CL3. Now q, II, B have the same expressions as above and 

/i(0) = c, /i(B/2) = a 2 , /J,(ni + B/2) = a\, /i(vri) = a 3 . 

The function fi(z) is real and finite only along the real axis, on which the differential has a double 
zero. As a result, the composition U(X) = /i(z(X)) varies between a 2 and c, with the same period 
as above and has cusps for U = c. See figure lb. 

The cases of other positions of c among the real roots (ai, a 2 , a 3 ) lead to either one of the above 
two types of solutions or to unbounded solutions. 



3 Periodic and solitary peakon solutions. 

In this section we describe different periodic and solitary solutions to the Dym type and shallow 
water SW equations obtained as limits of the traveling wave solutions and associated with various 
singularizations of the elliptic curves described in the previous section. In particular, we encounter 
periodic and solitary peakon solutions, i.e. solutions with discontinuous derivatives. 

In contrast to traveling wave solutions given above globally in a parametric form, peakon solu- 
tions can be given explicitly only in certain intervals. 



Periodic peakon solution for the Dym type equation. Consider again the integral ( p.l| ). 
Assuming k = —1, in the limit a 2 — > c the latter can be written in the equivalent differential form 

. dX = dX, or dX = dX', dX = (X-c)dX', (3.1) 

V«i — A (A — cjyai — A 
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X' being a new variable. Now putting A = a\ — v 2 , a± — c = a 2 and integrating the second 
differential in ([O]), we find that 

1, fu-a\ v , „, „, . l + e aX ' 

log )=X' + C, C = const, i.e., u = a 



a~~°\v + a) "l-e aX> ' 
Let us choose here C = iri/a. Then we have 

( e -aX'/2 _ e aX>/2\ 2 
X(X') - c = al-a\[ e _ aX , /2 + eaX , /2 J , or, equivalent^, 



X(X') - c = 4<9|, log t(X'), t(X') = e ~ aX '/ 2 + e aX '' 2 . (3.2) 
This function is 2-7ri periodic in aX and takes real finite values along the real axis. 



Next, using the expression (3.2), we integrate the relation dX = (X(X') — c)dX' and get the 
connection between X and X': 

e -aX'/2 _ e aX'/2 

X = 4d x , logr(X') + X = 2a e _ aX , /2 + eaX , /2 + ^o, = const. (3.3) 



Notice that as X' — ► — oo or +oo, X has finite limits which differ by 4a. Thus, expressions Q3.2| ) and 
(|3.3| ) give the solution U(X) = X(X'(X)) in a parametric form only in the interval [Xq— 2a, Xq+2o\. 
As follows from ( |3.2| ) and ( |3.1| ), at the endpoints of the interval, U = c and dU/dX = ±\/ a i — c - 

Now we define the global solution U(X) for all values of X by using periodic extension, i.e., by 
gluing an infinite number of pieces corresponding to Xq = AaN, N £ 7L in ( |3.3j ). At the endpoint 
of each interval the derivative of the solution changes sign resulting in a peak, and we obtain a 
periodic peakon solution. 



On the other hand, a direct integration of the first differential in (3.1) yields the following 
solution which holds between subsequent peaks 

X{X) = oi - X -{X - X ) 2 , X £[X - 2a, X + 2a]. 

Thus, in contrast to X(X') and X{X'), the profile of U(X) between peaks is not exponential, but 
a quadratic one. 

Periodic peakon solution for the SW equation. In a similar way, consider a limit of the 
periodic solution of the SW equation by putting 03 = c in (^S|). Then the curve £ becomes singular 
having a double point at fj, = c. Let £' be the corresponding regularized curve. The integral ( |2.9[ ) 
gives rise to the following differentials 

^ -dX, dfi =dX>, (3.4) 



y/(n- ai)(n - 02) {n - c)y/Jp - ai){/j, - a 2 ) 

where, as before, X = x—ct, dX = (fi—c)dX'. Integrating and inverting (ET1) we obtain respectively 

K x ) = 1 (e + e ) + o » ( 3 - 5 ) 
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^ ' ( e Z/2 + e - Z/2)2/( 0l _ c ) + ( e Z/2 _ e -Z/2)2/( a2 _ c ) ' (3,6) 

£ = V(c-ai)(c-a 2 )X'. 

The last expression is 27rz-periodic in Z and provides a parameterization of the regularized curve 
£'. Integrating (|3.6|) with respect to X', we find that 



v , (ai-a 2 )e z - (a x + a 2 ) + 2c - 2^ - ai)(c - a 2 ) . . 

A = log + const. (3.7) 

(ai - a 2 )e~ z - (oj + a 2 ) + 2c + 2-y/ (c - Oi)(c - a 2 ) 



Now suppose that 0,1,02, c are real and c < ai < a 2 or a\ < a 2 < c. As follows from (3.6), ( |3.7[ ). 
for real parameter X', the variables /U and X are real as well. As X' — * ±00, the function X(X') 
varies between different finite limits which we denote by X\ and X 2 , whereas, by ( |3.6| ), fi(X') — c 
tends to zero. As a result, the composition function U(X) = n(X'(X)) is defined only in the 
interval L\i,X 2 ], where it takes values in [c,a\] or [a 2 ,c]. According to (|3.4j) , at the endpoints of 



the interval, U = c, dU/dX = iy 7 (c — ai)(c — a 2 ). Like for the Dym equation, we define the global 
solution U(X) by periodic extension of /j,(X'(X)). This implies that U(X) has periodic peaks on 
the infinite interval (— 00, 00). 

The expression in (|3.5| ) provides a piece of the solution between two subsequent peaks. Solving 
the equation fi(X) = c or using ( |3.7| ) directly, we find the period of the peakon solution to be 



X 2 -X 1 = log 2c-(ai + a 2 ) + 2/(c-a 1 )(c-a 2 ) (3 g) 

2c - (ai + a 2 ) - 2-^/ (c - ai)(c - a 2 ) 

See figure lc. This completes the description of the periodic peakon solution. 

One can show that the other possible case: a\ < c < a 2 does not provide a real bounded 
solution. For details about algebraic geometric approach to describing n-peakon solutions see 
Alber and Miller [2000] and Alber et al. [2000a]. 



Solitons for the SW equation. Consider another possible degeneration of the integral (|2.9| 
assuming that a 2 = 03 = a and a, a\, c are distinct. Then we have 



(/i - c) ^ -dx, , it = dx\ 



(fj, - a)yj{n - ai)(n - c) (n - a)^J{ji - oi)(ju - c) 

X = x — ct + const, dX = (p — c)dX' . 



(3.9) 



The first differential has 2 pairs of simple poles on the rational curve 

{v 2 = (fj, — a\){(jL — c)} corresponding to \i = a and fi = 00 and a double zero for fj, = c. Therefore 
the inverse function U(X) again has branching of the form ( |2 . 2| ) . 

Integrating the second differential in (p.9|) and inverting, we obtain 



4 

a(X') -a = = = = , (3.10) 

(e z / 2 + e~ z / 2 ) 2 /(c-a) + (e z / 2 -e- z / 2 ) 2 /(a 1 -a) 



Z = \/ (a — a\)(a — c)X'. 
Then, integrating the relation between X and X' yields 



Y i \v'-li (ai - c)e z - {a x + c) + 2a - 2^(0 - ai)(o - c) 

A = (a — c)A + log = + const. (-J"- 1 ^) 

[a\ — c)e~ z — (ai + c) + 2a + 2^ (a — a\)(a — c) 

As before, let us consider possible real solutions assuming that a, a\,c are real. According to 
( 3.10| ), n{Z) is 27ri-periodic and it takes real values along the lines ImZ = Niri, JV£Z. 



10 



Case 1. a < a\ < c. The function /j,(Z) is real and finite only along the lines ImZ = iri + 2Nni, 
where it describes a smooth solitary wave tending to a as MX or MZ tends to ±00 and having the 
maximum a%. Along these lines the derivative dX/dX' = fj,— c is always real, negative and separated 
from zero. Therefore, the composition function U(X) = /i(A'(A)) gives a smooth solitary wave as 
well. 



Case 2. a < c < a\. Now p(Z) is real and finite only along the lines ImZ = 2Niri, where it again 
describes a smooth solitary wave tending to a as MX — > ±00 and having the maximum c. On the 
other hand, along these lines the derivative dX/dX' = fj, — c is always real and negative except 
when fi = c. At this point dX/dX' has a double zero. As a result, for fi = c the derivative dll/dX 
blows up and the composition function U(X) = fi(X'(X)) describes a solitary cusp (cuspon). 



Solitary peakon for the SW equation. As seen from ( p.8[ ), in the limit a± = a>2 = a the period 
of the above peakon solution tends to infinity. This gives us a solitary peakon solution, which can 
also be regarded as the separatrix between the smooth soliton and cuspon solutions. Indeed, in 
this case the algebraic curve £ has 2 double points with fi = a and /U = c. Its reauthorization £ 
consists of two disjoint copies of P = {/i} U 00. On them the differentials (|3.4j) take the following 
simple form 

d/i -IdX, - ^ - = ldX', (3.12) 



fx — a ' (n — c)(/i — a) 

where, as above, dX = (pi — c)dX' and / = ±1, regarding to the copy of P. First, we suppose I = 1. 
After integration in \i and inversion this yields 

p(X)=a + e x+c \ C 2 = const, = °~ C ^)x> • ( 3 - 13 ) 

The function fJ,(X') has period 2m /(a — c) and its fundamental domain gives a parameterization 
of the cylinder P \ {/i = a} \ {fi = c}. 

Integrating the differential relation between X and X' , we obtain 

f X ' (a - r)p- ( - a -^ X ' 

X ^ X ') = - / -1 (a C )x> dX ' = ~ Me (c -° )X - 1) + Ci, X' , d = const. 

Jx' l - e-y a - c ) x 

Suppose that a,c are real and, for definiteness, put a < c. The function n(X') is 2ni/(a — c) 
periodic, it has a pole along the real axis, whereas along the line ImA' = wi/(c — a) it varies in 
the interval (a, c). Let us put Cx = vri. Then X(X') is real along this line. We notice that as 
MX' — ► 00, A tends to —00 and /i tends to a, whereas for MX' — > —00, we have X — > and /i — > c. 
Thus the composition function n(X'(X)) is defined on the interval [— 00, 0] only, where it is also 
given by the first expression in ( |3.13[ ) for C2 = log(c — a). 

Now assuming I = —1, we obtain the same expressions with X replaced by —A. Let us choose 
the same values of the integration constants C\,C%. Then the composition function p(X'(X)) is 
defined on the interval [0, 00], where it has the form //(A) = a + (c — a)e~ x . As a result, the two 
copies of P give rise to two branches of the real continuous solution 

U(X) = a + (c — a)e~' x ', let, 

which has a peak at the origin. 
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Remark. As seen from the above peakon solutions, the amplitude of peaks coincides with the 
velocity c, which is a specific property of soliton propagation. 



4 Stationary quasi-periodic solutions. 

Stationary solutions provide profiles of the quasi-periodic wave solutions. For the sake of clarity, in 
this paper we restrict ourselves to the simplest nontrivial case g = 2. All the formulae and solutions 
below (except structure of real solutions) can be easily extended to the arbitrary (^-dimensional case 
which is only notationally more complicated. 



Stationary quasi-periodic solutions for the SW equation. According to the trace formula, 
in the genus 2 case we have 



5 

U(x,t) = m + fj, 2 - y^aj, (4.1) 



and equations (|1.3|) take the form 



H , , = at, 



n\dm + i4dfj, 2 = dx 



(4.2) 



2/Re(w) 
where 

Re(fj-) = —Kfj,(fj, — a\) ■ ■ ■ (fi— 05), ai, . . . , 05 = const. 

Here we suppose that all the roots of Re(fjL) are distinct. The variables fJ-i,fJ> 2 must be regarded as 
coordinates of points Pi = (ni,wi), P2 = (^2, w 2 ) on the genus 2 hyperelliptic curve T = {w 2 = 
Reifi)}- Equations ( |4.2j ) involve one holomorphic differential and one meromorphic differential of 
the third kind having a pair of simple poles at the infinite points oo_, oo + on T. Integrating ( f4.2[) , 
we obtain the mapping of the symmetric product to C 2 = (t, x) 

Pl [idfi ^ l' P2 ndfi^ ^ 

(4.3) 



P 2 v / i^) Jp 2y/R^L) 

+ / , , = X. 



Po 2^R 6 (fi) Jp 2V Rq{^) 

where Po is a fixed basepoint of the mapping. Notice that for Pi or P2 = oo_, oo + , we have x = 00. 
Let us fix a canonical basis of cycles A\ , A2 , Pi , B2 on V in a standard way (see, for example, 
Mumford [1978]). The mapping has four independent (over the reals) 2-dimensional vectors of 
periods of the above differentials along the cycles. In addition, it has one extra period vector 
corresponding to a homology zero cycle around oo_ or oo + . As a result, the mapping ( [4,3j ) has 
5 period vectors in C 2 hence its inversion is not well defined: there do not exist meromorphic 
functions on C 2 with five periods. In particular, U(x, t) is not a meromorphic or single valued 
complex function of t, x. 
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In order to describe properties of U(x,t), we, first, fix time by putting t = to (dt = 0) and 
consider stationary solutions U(x,to). Introduce a new coordinate x' such that 

dx = Hi/J.2 dx' . (4-4) 

Then equations ( |4.2| ) lead to the Abel-Jacobi mapping of to the Jacobian variety Jac(r) of T, 
which includes holomorphic differentials only: 



Po 2 v / i^(70 Jp 2^Re(ji) 

Pl fidfi , f P2 ndfi _ (4.5) 

ui = x + const, ii2 = const, 

where u\,U2 are coordinates on the universal covering C 2 of Jac(r). 

Let uj\ , CJ2 be the dual basis of normalized holomorphic differentials on V with respect to the 
above choice of cycles and z\,z 2 be the corresponding coordinates on the universal covering of 
Jac(r): 

d u + dun _ d 2 i + d 22 ^ , 

uji = . du, uj 2 = , du, 

2^BM (46) 

Z\ = d\\U\ + d\ 2 U 2 , Z2 = d2\U\ + ^22^2- 

Here the normalizing constants d are uniquely determined by the conditions § A . uJj = <5y . 

Recall that the standard theta-function related to a Riemann surface of genus g and theta- 
functions with characteristics a = («i, . . . , a g ), (3 = {(3\, . . . , (3 g ) £ W g have the form 

9(z\B)= expQ(BM,M) + (M,z)V 

M6Z9 ^ ' 

9 9 

(M, z) = Y^ MiZi, (BM, M) = BijMiMj, (4 ?) 



a 
P 



i=l i,j=i 
(z\B) = exp{(Ba, a)/2 + (z + 2m/3, a)}6(z + 2iri/3 + Ba|B), 



B being the g x g period matrix of T. In the sequel we shall omit it in the notation. 

Now we choose the basepoint Po of the mapping ( |4.5[ ) to be the last Weierstrass point (as,0) 
on r. Then, according to the trace formula for even order hyperelliptic curves (see e.g., Clebsch 
and Gordan [1866], Dubrovin [1981]) 

U = Mi + y«2 - J2 a 3 = const ~ d w log , (4- 8 ) 

0[5\(z + q/2) 

T f°° + 
J CO — 
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where, in view of the normalizing change ( |4.6|) , Z\ = d\\x' + const, zi = di\x' + const, dw is the 
derivative along a tangent vector W of T C Jac(r) at oo+, namely, in the coordinates (111,112), W = 
(0, 1) T , and in the coordinates (z x , z 2 ), W = (d 12 , d 22 ) T - Finally, 5 = {6", 5') T , 5", 5' E ^/Z^ is the 
half-integer theta-characteristic corresponding to the vector of Riemann constants (see Mumford 
[1983]). For the chosen standard canonical basis of cycles and the basepoint Pq = (05,0) 

5' = (1/2, . . . , 1/2) T , 5" = (g/2, (g - l)/2, . . . , 1, 1/2) T (mod 1). (4.9) 

Thus, in our case 

1/2 
1/2 1/2 

The function U(zi,Z2) is meromorphic on Jac(r) and it has simple poles along 2 translates of the 
theta-divisor 9 = {6(z) = 0} C Jac(r): 

©_ = {9[5]{z - q/2) = 0}, G + = {6[5}(z + q/2) = 0}, 

which are tangent to each other at the origin {z = 0}. Thus, U(zi(x'), Z2{x')) is a quasi-periodic 
function of the complex variable x'. Notice that a quasi-periodic genus 2 solution of the nonlinear 
mKdV equation has the same form. 

We also notice that the point Eq = (// = 0,w = 0) is a Weierstrass (branch) point on T. 
Then, following Clebsch and Gordan [1866], we have the following expression for the symmetric 
polynomial 

2 [S + Vo](z) 
6[8}(z + q/2)0[6]{z - q/2) 

where t]q is the half- integer theta-characteristic corresponding to the branch point Eq: 

rE 

Vo = (Vo,Vof> / (^1,^ = 27^ + 13^ eC 2 . (4.11) 
JPo 

Thus, the product is a meromorphic function on Jac(r) having simple poles along 0_, Q + 
and a double zero along another translate of the theta-divisor 0, ©o = {9[5 + T]q](z) = 0}, passing 
through the origin and intersecting each of the translates 0_, 0+ at two points. The translate ©o 
can also be interpreted as the image of the curve T itself under the Abel-Jacobi mapping (| 



P rE 
(U!1,U) 2 ) T + / 



©0 = { I (UI,IV2) T + / (UJI,UJ2) T 



Per 



It follows from (4.4), ( 4.1Q ) that generically the derivative of the function x(x') is equal to and 



that it has a double zero each time when the complex x'-flow intersects ©0, i.e., when 6 [6 + t]o]( z ) 
vanishes, except the points where the flow is tangent to ©0, i.e., when 6[S + t]q](z) has a higher 
vanishing order in x' and \1\H2 too. This takes place only at the origin of Jac(r). Since at the 



origin the solution (4.8) blows up, we conclude that for bounded solutions the function \i\\i2 may 
have only a double zero and x(x') a simple zero in x' . 

On the other hand, in view of the second equation of ( |4.3j ), the original variable x is a sum of 
Abelian integrals of third kind. Introduce the normalized differentials of third kind ^oo_oo + on T 
having poles at oo_,oo + with residues ±1: 

^oo_oo + = 7== +h 1 u 1 + h 2 2 , (4.12) 
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where hi, h% are normalizing constants specified by the conditions for 0, 00 _ OQ+ to have zero ^4-periods 

on r. 

According to Clebsch and Gordan [1866], 

Pl f P2 6\8\(z 4- a/2) 

^oo_oo + + / ^oo_oo + = log a \J.) ^ qi > + const. (4.13) 

Then, in view of the second equation in ( |4.2| ) and ( [4.12D , we get 



, „ , 6[5](z + q 2) , 

x(x ) = log — — — - hizi - h 2 z 2 + const, 

6[5\(z - q/2) 

%l = dux + const, Z2 = d,2ix' + const. 



(4.14) 



As a result, we expressed the stationary quasi-periodic solution U and the argument x in terms 
of the auxiliary complex variable x' . The algebraic geometrical structure of the general solution 
U(x, t) and the behaviour of real solutions will be considered in the next sections. 

Stationary quasi-periodic solutions for the Dym equation. Now we pass to the Dym 



equation (lL2j and seek its solutions again in the form (4.1). In this case the variables [ii,[J-2 again 



change according to equations of the form (4.2) with the only difference being that the order of the 



polynomial defining the corresponding hyperelliptic curve is odd: 



2^/RM 

lAdyi (4^2 ^ dx (4.15) 

2yTR^j 2^fR^2) X ' 

-R5(m) = -nfi{fj, -ai)---((j,- a 4 ), 

hence the corresponding hyperelliptic curve T = {w 2 = R^,(fi)} has just one infinite point oo. As a 
consequence, the equations ( [4.15| ) contain one holomorphic differential and one differential of the 
second kind. 

As before, we first consider stationary solutions by putting t = to (dt = 0) and assuming k = 1. 



Notice that under these conditions, (4.15) has the same structure as quadratures for the Jacobi 
problem on geodesies on a triaxial ellipsoid Q, where fii,(J>2 play the role of ellipsoidal coordinates 
on Q, parameters oi, 02, 03 the squares of the semi-axes of Q, 04 the constant of motion, and x the 
length of a geodesic. 



Under the change of parameter (4.4), we arrive at the Abel-Jacobi mapping 



Pl dfi r P2 dp _ 

P 2 x /iU^) + 7p 2^/R^l) ~ UU 

Pl fidfi f P2 udu (4.16) 

+ / — : : = U2, 



P 2^/EM J Po 2y/R 5 fa) 
ui = x + const, u 2 = const 



15 



This change was first made by Weierstrass [1878] in order to find the theta-functional solution for the 
geodesic problem (see also Cewen [1990]). Next we introduce normalized holomorphic differentials 
a>i,o>2 on r and coordinates Z\,Z2 on the universal covering of Jac(r) according to ( |4.6| ) and, in 
addition, the normalized differential of the second kind having a double pole at oo 

nW = -ife + /^ 1 + /^ 2 . (4.17) 
2 v / i? 5 (/i i ) 

Similarly to (4.12), the constants h[, h' 2 are uniquely defined by the condition that Q^a have zero 
^4-periods on T. 

Then, instead of the expressions ( |4.8j ), and ( [4.10 ), we have (see, e.g., Dubrovin [1981], Dubrovin 
et al. [1985]) 

U(x') = m +/x 2 = const - d v 9[5](z), 

(4.18) 

z\ = d\\x' + const, Z2 = d2\x' + const, 

where dy is the derivative along the tangent vector V of T GJac(r) at oo: V = (d\2, c?22 ) T 7 and, 
respectively, 

9 2 [5 + m }(z) , A , 

= « — fl n m/ ; — , « = const, (4.19) 
6^[<5J(z) 

where the characteristic 770 is defined in ( (4.11| ). In addition, in contrast to ( P^D , the sum of Abelian 
integrals of second kind has the form 

[ Pl fi&> + r 0« = const - 9 y log 6[5](z). (4.20) 

Comparing this with ( |4.17| ) , we find that the analog of the relation ( 4.14| ) between the parameters 
x and x' has the form 

x(x') = — dy log 0[<5](z) — h\Z\ — h,2Z2 + const, 

(4.21) 

z \ = d\\x' + const, z 2 = d 2 \x' + const. 

This expression can be regarded as a 2-dimensional generalization of the Weierstrass zeta-function 
in Q. 

Thus, we have expressed the stationary solution U and the argument x in terms of the auxiliary 
complex variable x'. Various types of real solutions defined by the above expressions will be 
considered in Section 6. 



5 Time-dependent quasi-periodic solutions. 

The solutions for the SW equation. In order to obtain general time-dependent solutions 
U(x,t) of the SW equation given by (4.1), one has to invert the mapping (f4.3p. However, as 



already mentioned, the problem of inversion is unsolvable in terms of meromorphic functions. 

To describe the structure of general solutions, let us first consider a divisor of 3 points Pj 
(/Xi, Wi), i = 1, 2, 3 on T \ {oo_, oo + } and the following extended equations 



3 3 

dm _ An ^ Hidfii _ ^ + ^ fx'fdfii 



dy, y = dt, y =dx, (5.1) 

-2/R^) f^2^RM ^2^RM 
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including the extra variable y, two holomorphic differentials and one differential of the third kind 



on r. The latter are linear combinations of the normalized differentials 0)1,0)2, ^±00 defined in (4,6) 
and ( 4.12j ), According to Clebsch and Gordan [1866], Fedorov [1999], equations ( |5,1|) describe a 
differential of a well defined mapping of the symmetric product (T \ {oo_, oo + })( 3 ) to generalized 
Jacobian variety Jac(T, oo±), a noncompact algebraic group represented as the quotient of C 3 by a 
lattice A generated by five vectors of periods of the differentials 0)1,0)2, ^±00 on T. Topologically, 
Jac(T, oo±) is the product of the 2 dimensional variety Jac(T) and the cylinder C* = C \ {0}. An 
analytical and algebraic-geometrical description of generalized Jacobians can be found in Clebsch 
and Gordan [1866], Belokolos et all [1994], Fedorov [1999], Gavrilov [1999]. 

Let Zi, 02 j Z be coordinates on the universal covering of Jac(T, oo±) such that 

3 ,p t 3 , Pt 3 „ Pt 

} ) I V! = Zl, / u 2 = z 2 , } y I ^±oo = Z, (5.2) 



where, as above, Pq = (05,0). Then, according to (|4.6| ), ( [4.12D , 

Z\ = dny + d\ 2 t + const, z 2 = d 2 \y + d 22 t + const, , , 

Z = x + hi(d\iy + di 2 t) + h 2 (d 21 y + d 22 t) + const. 

The problem of inversion of Abel-Jacobi mappings including differentials of the third and second 
kind is solved in terms of generalized theta-functions which are finite sums of products of customary 
theta-functions, rational functions, and exponentials (see Ercolani [1989], Fedorov [1999], Gagnon 
et al. [1987]). To invert the mapping (|5.2| ), we shall make use of the following theta-functions 



9(z, Z) = e Z l 2 9{z + q/2) + e- Z ' 2 6{z - q/2), 

(5.4) 

9[r,](z, Z) = e z / 2 9[ v ](z + q/2) + e~ 2 '9[ V ](z - q/2), 

z = (z 1 ,z 2 ), q=(qi,q 2 ) , qi = td 1 , q 2 = / tu 2 , 

Joo— Joo_ 

where 9(z), 9[r]](z) are customary theta-functions associated with the curve T with half-integer 
theta-characteristics 77. Like 9[rj](z), generalized theta-functions have a quasi-periodic property: a 
shift of the argument (z, Z) by any period vector of the generalized Jacobian results in multiplication 
of #[r/](,z, Z) by a constant factor. 

Now consider the dissection r of T along the canonical cycles A± , A 2 , B\ , B 2 , which is a one- 
connected domain having the form of an octagon. In addition, we cut V along the paths joining a 
point O on the boundary dT of V to the points oo_, 00+. On the obtained domain V we introduce 
the sing le-valued function F(P) = 9[5}(A(P) - (z, Z) T ), where 

A(P) = ( [ Ol, f LD 2 , [ , 

\Jp JPo JPo / 



and the characteristic 5 is defined in (4.9). Then the following analog of the Riemann theorem 



holds (see e.g., Fedorov [1999], Gagnon et al. [1985]). 

Theorem 5.1 Let the coordinates z,Z be such that the function F(P) does not vanish identically 
on T' . Then it has precisely 3 zeros Pi, P 2 , P3 giving a unique solution to the inversion of the 
generalized mapping ([57 
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Now let us consider the logarithmic differential n(P)d log F(P). By Theorem 5.1, the sum of the 
residues of its poles in the domain V' equals /*(Pl) + ^{Pz) + M(-Fb)- Applying the residue theorem, 
after calculations, we get the following compact "trace formula" 

e z 6[5\(z + q)+e- z 9[5\(z-q) 
Ml + M2 + /*3 = const ^ 6[5}(z) ( ) 

with the characteristic 5 specified in fl3~9|). 

The principal difference between the extended mappings ( |5.1| ) or (^^) and the system ( [4.3|) is 
that the latter contains only 2 points on T \ {oo_, oo + }. On the other hand, ( |5.1| ) reduces to ( |1.2| ) 
by fixing P3 = Pq (/Z3 = 05, dfis = 0). Under this condition, ( |5.2|) describes the embedding of 
the symmetric product (r \ {oo_, oo + })( 2 ) into Jac(r, oo±). Its image is a 2-dimensional nonlinear 
analytic subvariety (stratum) W 2 . Like the generalized Jacobian itself, it is a noncompact variety. 



Remark 5.1. In case of customary Jacobian varieties, the corresponding nonlinear subvarieties 
and their stratification have been studied in Gunning [1972] and Vanhaecke [1995]. Such varieties or 
their open subsets often appear as (coverings of) complex invariant manifolds of finite-dimensional 
integrable systems (see Vanhaecke [1995], Abenda and Fedorov [1999]). 

It follows from the above that on the stratum W% the variables zx,Z2,Z play the role of excessive 
(abundant) coordinates, hence they cannot be independent there. The analytic structure of W 2 is 
explicitly described by the following theorem (see e.g., Fedorov [1999], Gagnon et al. [1987]). 

Theorem 5.2 The subvariety W 2 C Jac(r, 00 ±) coincides with the zero locus of the generalized 
theta-function: 

W 2 = {e Z/2 6[5](z + q/2) - e~ z / 2 9[5](z - q/2) = 0}. (5.6) 



On the other hand, in view of relations (|5.3| ), the coordinates z,Z are linear functions of the 
variables x, t, and y. Thus, equation ([O]) can also be regarded as a constraint on them. It follows 
that on fixing P3 = Pq, y becomes a transcendent function of x,t. 

Now we notice that the sum fj,i+ /X2+Q5 = fJ-(Pi) + M-P2) + fJ>(Po) coincides with the restriction of 
the total sum //(Pi) + //(P2) + /*(^ > 3) ) as a function on Jac(r, oo±), onto W%. Then, using expression 
(|5.5|), we conclude that the 2-phase solution of the SW equation has the form 

rr, e z 6[5](z + q) + e- z 6[5](z-q) 

U(x, t) = const — q > 1 n ^ (5.7) 

9[5\{z) 

z\ = duy + d 12 t, z 2 = d 2 iy + d 22 t, Z = x + hi(dny + di 2 t) + h 2 (d 21 y + d 22 t), 

where the extra variable y depends on x, t according to (5.6). As a result, we arrive at the following 
algebro-geometric description of motion: 



The x-flow (t-flow) defined by equations ( j.i ) evolves on the nonlinear variety W 2 C Jac(r,oo-i-) 
in such a way that y is a nonlinear transcendent function of x (respectively, oft). In this sense the 
flow is nonlinear. 

We emphasize that the solution U(x, t) is neither meromorphic in x, nor in t. 
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Remark 5.2. Let us consider the x-flow by putting t =const. It turns out that, up to an additive 
constant, the extra variable y can now be identified with the auxiliary variable x' introduced in 
( ]4.4j ) when we considered stationary solutions. Indeed, in view of ( |5.3D , in this case the condition 
in ( |5.6D becomes 

Z = x + h\Z\ + h 2 z 2 + const = log 7L-L7— — + const, (5.8) 

0[<$](z + 9/2) 



which is equivalent to the relation ( |4.14| ) between x and x'. In view of (5^8) and the addition 



theorem for theta- functions, the solution ( |5.7| ) reduces to the stationary solution 

In contrast to x, the parameter t enters both expressions for z and Z in Q5.7| ). Therefore, in the 
case of the i-ffow, t cannot be explicitly expressed in terms of y as in the case of the x-Qow. This 
implies that solutions U(xQ,t), xq =const must have different properties in comparison with (|4.8[). 



Remark 5.3. We notice that the subvariety W 2 of the generalized Jacobian becomes linear in 
rare cases when the curve T enjoys some nontrivial involutions, i.e., when it can be regarded as a 
covering of an elliptic curve. (Various examples of the involutions can be found in Belokolos et al 
[1994].) In such cases U(x,t) becomes a meromorphic function of its arguments. 

The solutions for the Dym equation. Now we proceed to the problem of inversion of the 
reduction ( [4.15|) of the Dym equation which is related to the odd order hyperelliptic curve T = 
{w 2 = i?5(^)}. As in the case of the reduction of the SW equation, in order to describe the function 
U(x,t) = n\ + /U2, we first consider an "excessive" divisor of 3 points Pj = (fii,Wi), i = 1,2, 3 on 
r \ {00} and the extended equations 

^2^R^h) ^2/R^) ^2^RM (5.9) 

Rbin) = -kh(h - a\) ■ ■ ■ (/i - a 4 ), 

including 2 holomorphic differentials and one differential of the second kind having a double pole 
at 00 G r. They are linear combinations of the normalized differentials o>i, office defined in (4.6) 
and KT%. 



In contrast to (5.1), equations (|5.S|) describe a differential of a well defined mapping of the 



symmetric product (F \ {oo})( 3 ) to the generalized Jacobian variety Jac(r, 00), the quotient of C 3 
by the lattice generated by four period vectors of the differentials u>i,u>2, on T. Topologically, 
this variety is a product of the 2 dimensional variety Jac(r) and the complex plane C (see Clebsch 
[1866], Gavrilov [1999]). 

Let us introduce coordinates Z\, Z2, Z by the mapping 

3 ,.p 3 ,.p 3 ,.p 

y r^ = Zl , y jco 2 = z 2 , yrnw = z (5.10) 

i=l Je o i=i JEo i=1 JE 

with the basepoint Eq = (0, 0) (we cannot choose the basepoint to be 00 as in the previous section, 
since it is the pole of )• Next, comparing (4.6), ( 4.17 ) with (5.9), we find the following relations 



z\ = d\xy + d\it + const, z 2 = d 2 \y + d 22 t + const , . 

Z = x + h[(duy + di 2 t) + h' 2 (d 2 iy + d 22 t) + const. 
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Like Q5.2D , the mapping (5.10) is invertible in terms of meromorphic functions. The inversion 
problem is solved by means of the following rational degeneration of the customary theta-function 

e(z,z) = ze[5](z) + d v o[5}(z), (5.12) 

where dy is defined in ( |4.1&j ) (compare with the generalized theta- functions Q5.4| )). Like ( |5.4[ ), the 
function ( 5.12|) enjoys the quasi-periodic property. 

Consider again the dissection r of V and cut it along a path joining a point O on the boundary 
dT to 00. In the obtained domain we introduce the single- valued function 

F(P) = (z- / P ft«W] (z- f°°Q- [ P u)+d v e[5] (z- f°u>- f q). 

V JEo J \ JE JEo / \ JEo Je J 

Using a modification of Theorem 5A and calculating the logarithmic differential n(P) dlogF(P), 
we obtain 

Ml + M2 + A«3 = const - {Z + dy log 9{8 + Vo}( z )) 2 ~ d v log 9[8 + i] ](z) 

2 2Zd v e[6 + m ](z)-e^e[6 + i 1 o](z) (5 - 13) 

— const — Zj -77 t- — r , 

9[8 + m ]{z) 

where 770 = (r]Q,r]' ) T 6 |Z 2 /Z 2 , such that 2nir)Q + Bt/q = W2) T - Now, similarly to the case 

of the SW equation, we fix P3 = Eq = 0, ^3 = 0) in the mapping ( 5.10Q . In this case its image 
becomes a 2-dimensional nonlinear noncompact analytic subvariety W2 C Jac(T,oo). Comparing 
the third sum in (530) and expression (4.20), we find 

W 2 = {Z + const + dv log 9[S + r] ](z) =0}. (5.14) 

Finally, taking into account the trace formula ( |5.13| ) , we conclude that the solution of the Dym 
equation has the form 

U(x,t) = = const - d v log 9 [8 + rj Q ](z), 

(5.15) 

z\ = d\\y + d\2t + const, z 2 = c?2iy + ^22^ + const, 

where the extra variable y depends on x,t in a transcendental way according to the constraint 
(|5.14j ) and the expression for Z in (|5.11| ). The solution U(x,t) is not meromorphic with respect to 
its arguments. 



Remark 5.4. As in the case of SW equation, the stationary solutions for the Dym equation given 
in the previous section can be obtained from (5.15) by putting t =const. Then y can be identified 
with the auxiliary variable x 1 defined in ( iA ) and the condition in fl5,14| ) becomes equivalent to the 
relation (4.21|) between x and x' . As a result, (|5.15) gives precisely the stationary solution (4.18). 



6 Real bounded stationary 2-phase solutions. 

In this section we impose reality conditions on the stationary complex solutions obtained in Section 
4. 

Let a be the antiholomorphic involution on a hyperelliptic curve T = {w 2 = P(p-)} of genus 
g. The part of T which is invariant with respect to a is called the real part r(M). On the plane 
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M 2 = (5fyx, "Siw) it is either the empty set or a union of ovals. By the Abel-Jacobi mapping, the 
involution a lifts to Jac(T). By Jacu(T) we denote the real part of Jac(r) that is invariant under a. 
One can show that the elementary symmetric functions of the variables fix,...,fj, g take real values 
on JacR(r) and only there. 

Theorem 6.1 (Comessatti [1924])- Let s be the number of connected components o/T(M) and L 
be the number of connected components of JacR(r). If s ^ 0, then L = 2 S ~ 1 . If s = 0, then L = 1 
provided the degree of R(fJ-) is even and L = 2 in case the degree is odd. 

Shallow water equation. Suppose all the roots of the polynomial Re(n) in fl4.2|) arising in the 
reduction of the SW equation are real, i.e., T(M) consists of 3 ovals about the segments [0, ei], 
[e2, es], and [e^, es]. By Theorem 5.1, JacR(T) has 4 connected components. They are characterized 
by the following behavior of ^-variables, which reflects in different properties of real stationary 
solutions U(x,to): 

Case 1. The variables H\,H2 are real and fi\ G [02,03], M2 £ [04,05]- The sum U = fi\ + //2 
is thus a real quasi-periodic function of x' having no poles and zeros. The product //1//2 has the 
same properties. Geometrically this means that the corresponding component of JacR(r) does not 
intersect the translates 0_,O+, and Bo- In view of ( f4.4p , x(x') and x'(x) are strictly monotonic 
real functions. Therefore the composition U(x,to) = U(x'(x)) is a quasi-periodic regular function. 

Cases 2,3. fj,% G [0,ai], whereas //2 £ [02,03] or [04,05]. The function U(x') has the same 
properties as above, whereas /X1//2 does not blow up, but has zeros. As found in Section 4, the 
derivative dx/dx 1 vanishes with a second order with respect to x' as one of the /x- variables vanishes. 
This happens when the real x'-fiow on the considered components of JacR(T) intersects 0o(K) = 
Oo H Jacx(r). If follows that at this moment the derivative dU/dx = dU/dx' ■ dx' jdx blows up and 
the graph of the function U(x, to) has an inflection point with vertical tendency line. 

In addition, when (^1,^2) = (0,02) or (0,04), i.e., when the real x'-flow passes a half-period on 
JacR(r), the function U(x') has an extremum, which implies that the graph of U(x,to) has a cusp. 
Due to quasi-periodicity of the flow, U(x,to) has an infinite quasi-periodic sequence of cusps. 

Case 4. Now the variables fii,^2 are complex conjugated. Using equations fl4.2| ), we show that 
they cannot reach real axis. It follows that the product /X1/U2 is always nonzero and U(x,to) is again 
a quasi-periodic regular function. 

In a similar way one can show that when some of the roots of R^ifi) are complex conjugate, the 
qualitative behavior of the real solution U(x, to) coincides with one of the above four cases. 

In a forthcoming paper we will consider different singular limits of the quasi-periodic solutions 
when the spectral curve becomes singular and its arithmetic genus drops to zero. The solutions 
are then expressed in terms of purely exponential tau-functions and, in the real bounded case, they 
describe an interaction of the two smooth solitons or cuspons, or a quasi-periodic train of peakons 
tending to a periodic one at infinity. 
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7 Peakon-soliton solutions and elliptic billiards. 

In this section we continue studying degenerate solutions of the Dym equation. Consider another 
possible confluence of roots of the polynomial -Rs(p) in ( 4.15| ) by putting 

ai = 0, a2 = 03 = 6, 04 = a, k = 1. (7-1) 

As before, we first analyze stationary solutions (x-flow) by setting dt = 0. Passing to the new 
variable x' according to the change 



dx = P1P2 dx', (7-2) 



from flOD 

we get 



dpi , dp 2 



+ / z = -feds', 



Hl\/a- pi /U 2 V« - P2 

dpi dp 2 = Q 



(7.3) 



(pi - 6)Va - Pi (P2 - &)V a - P2 



Let us introduce the following normalized differentials of third kind on the Riemann surface P = 
= a — p} that have pairs of simple poles Q± , Q± and Q% , Q 2 with p = and p = 6 respectively 

gidp ^dp , - 

p£ (P - &)£ 



Then equations ( |7.3|) give rise to the generalized Abel-Jacobi equations 

fii+ / fii = *i, / ^2+/ ^2 = 22, Pi = (jH,ti), (7.4) 
P ^Po M JPo 

z\ = -b(3ix' + zio, «io, 22 = const, (7.5) 

where we put Pq = (a,0). These describe a well defined mapping of the symmetric product 
(P \ {Qi,Qf, Q 2 ,Qi}) {2) to the generalized Jacobian Jac(P, Qf, Qf). 

As a result of inversion of ( [7.4|) , one finds the following expressions for symmetric polynomials 
of pi and P2 

Pi +P2 = U(z 1 ,Z 2 ) = d^r log t(z 1 ,Z 2 ) +01 ~Pi 

1 ( e ^/2 + g-2 2 /2\2 

P1P2 = -TZ-d zl d w logT(z 1 ,z 2 ) or pip 2 = iab T , 

bpi t 2 (z 1 ,z 2 ) 

d d 

dw = ^f3i- \-2f3 2 ——, z\ = -bflix + z w , z w , z 2 = const, 

az\ oz 2 

and t(zi,z 2 ) is the 2-dimensional tau-function with 01,0:2 replaced by (3\,(3 2 . The latter admits 
decomposition in the following sum of one-dimensional tau-functions 

r( Zl , z 2 ) = e z ^ 2 r(z 2 + q/2) - e~ z ^ 2 T{z 2 - q/2), 

t(z 2 ) = e 2/ +e 21 , q = log 

V Pi + P2 



where 
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Lastly, by using the second expression in ( |7.6[ ) and the relation between z\ and x' in (7.5), we find 

x = J /xi/i2 dx' = &w log r (z\, z%) + const 
or, in view of decomposition ([7,7|), 

/ a = fl e 2l/2 r(^2 + g/2) + e~^/ 2 r(z 2 - g/2) 
x ^ ) n eZl /2 T ( Z2 + q / 2 ) _ e-*i/ 2 T(z 2 - q/2) 
e^' 2 d Z2 r{z 2 + q/2) - e-^/ 2 a z2 r(z 2 - g/2) 
+/?2 e^/ 2 r(z 2 + g /2)-e-W 2r (, 2 - g / 2) + ^ (7 " 8) 

z\ = —b(3\x' + zio, Z2 = const. 



Remark 8.1. As mentioned in Remark 4.1, equations ( 4,15| ) with dt = describing the quasiperi- 



odic stationary solutions have the same structure as quadratures for the geodesic motion on an 
triaxial ellipsoid £ (more generally, a quadric) in M 3 . Parameter x plays the role of length of a 



geodesic. Under the limit (7.1) one of the semiaxes of £ tends to zero whereas the geodesic motion 



passes to the asymptotic billiard motion inside the ellipse 

£ = {Xf/a + Xl/b = l} Cl 2 = (Ii,I 2 ). 

Geodesies on £ transform to straight line segments passing through a focus of the ellipse between 
each subsequent elastic reflections (impacts) along £. As x — ► ±oo, the billiard motion tends to 
oscillations along the bigger axis of the ellipse. Now the variables /J>i,fi2 play the role of elliptic 
coordinates in M 2 such that 

2 _ (q-/ii)(a-/i 2 ) y2 _ (b - m)(b - fi 2 ) 

A l — 1 > A 2 — 7 • 

a — b b — a 

Along £ one of the variables equals zero. 

It follows that equations ([7l|) can be regarded as describing the straight line motion of a point 
mass inside £. When the point meets the ellipse, one of the /u- variables, say hi, vanishes, and the 
corresponding point Pi = , \J R(fi\)) on the Riemann surface P coincides with one of the poles 



Qi,Qi of the differential Q\. Then, as follows from the mapping (7A), z\ and x' become infinite. 

On the other hand, as x' ',zi — ► ±oo, the second expression in ( |7.6| ) vanishes, whereas the first 
one has finite limits giving the values of fi 2 at the subsequent impact points. The variable z 2 plays 
the role of a constant phase defining position of the segment between the points. 

According to (|7^), (TFj 



(7.9) 



and 



I \ a i a d ^ T ( z 2 ~ q/2) 

t(z2 - q/2) 

x(oo, z 2 ) = Pi+ 02 9z2 J Z2 , + !L 2 + c x , c x = const. 
t{z 2 + q/2) 



u( , B 2 R 2 

U{-OQ,Z 2 ) - -z _ ro — - „ + ft 



(/?i - /3 2 ) 2 e-^ + 09i + /? 2 ) 2 e^ + 2(/? 2 - /f ) 
^ ( °°' Z2) - (/3i-/? 2 ) 2 e- + (/3 1 +/3 2 ) 2 e--+2(/3 1 2 -/? 2 ) + A ~ ^ 2 " 



(7.10) 



23 



Notice that x(oo, z 2 ) and U(oo, z 2 ) have the same value as x{— 00, Z2+q , )+2/3i and U(— 00, 2:2+9) 
respectively. All this results in the following algebro-geometrical description: As the point mass 
inside £ moves from one impact to the next one, the point P\ on P moves from the pole to Qf. 
At the moment of impact, Pi jumps from Q± back to Q^, whereas the phase Z2 in ( [7.6| ) increases 
by q. Then the story repeats. 

Using this property, by induction, from ( 7.1C| ) the elliptic coordinates of the whole sequence of 
impact points are found in form 

„ =0 u = wm-pp 

^ ' ^ 2 (/3i-/3 2 ) 2 e^ + (/3i+/32) 2 e-^+2(/? 1 2 -/?| ) P*> (7.11) 

Z2N = Z20 + Nq, 

N £ N being the number of impact and the constant Z20 is the same for all the segments of the 
billiard trajectory. 

In addition, from ( fi 7 .^ ) we find the length of the iV-th segment of the billiard trajectory in form 

e <?/2 _ e -q/2 

x{oo, z 2N ) ~ x(-oo, z 2N ) = 2/3i + 2/3 2 r /2 -73 , (7.12) 

exp(z2Ar) - exp(-z 2 iv) + eil z - e II* 



Z20 being the same as in (7.11). 



According to the trace formula U (x, to) = /ii + /42, expressions ( [7.6|) , ( ffTq) provide us stationary 
peakon solutions of Dym equation in a parametric form. Here the phase Z2 must be regarded 
as a certain function of to- Namely, the expressions describe one piece of the profile U(x,to) 
corresponding to trajectory of the point mass between subsequent impacts. The other pieces are 
obtained by changing the phase z 2 in (|7.6|) , (frig) by q and adding 2/3\ to x. The pieces are glued 
at peak points, where the spatial derivative of U changes sign and which correspond to impacts in 
the billiard problem. The profile U(x, to) thus consists of an infinite sequence of peaks and knots 
between them. For this reason we call this solution the soliton-peakon solution. Notice that, in 
contrast to exponentials profiles of the functions U(x'), x(x'), any piece of U(x,to) has quadratic 
profile, as will be explained below. 

The heights Un of N-th. peak is given by ( [7. lip. T he distance between subsequent peaks along 
x-axis is a quasiperiodic function of given by ( 7,12; ). 



We emphasize that, in contrast to the peakon traveling wave solution considered in Section 3, 
now the rc-distance between subsequent peaks is not constant. However, as seen from (0, for 
N — > ±00 it tends to 2/3i, the doubled bigger semi-axis of the ellipse, whereas the pieces tend to 
identical ones corresponding to periodic billiard motion along Ai-axis. 



Remark 8.2. Expressions (7.6), (ff^) describe an asymptotic motion of an elliptic as well as a 
hyperbolic billiard. In the first case the initial phase z\o is pure imaginary whereas in the second 
case it is real. According to the trace formula, the hyperbolic billiard corresponds to unbounded 
stationary solutions of HD equation, which is out of interest of this paper. 



Remark 8.3 . The above considerations can be extended to multi-dimensional case. Namely, 
following similar approach one can consider genus g solution of HD equation described by equations 
( ]1.3| ), then its asymptotic stationary limit which is related to a generalized Abel-Jacobi mapping 
including g meromorphic differentials of 3rd kind. Consider a billiard inside a (/-dimensional ellip- 
soid. Then such a limit solution corresponds to asymptotic billiard trajectories that intersect g — 1 
focal quadrics of the ellipsoid between any subsequent impacts. The resulting stationary solution 
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U(x,to) consists of an infinite series of peaks and between each subsequent peaks there are g — 1 
knots. 

In order to study time-dependent soliton-peakon solutions, we consider the system ( |4 . 1 5| ) under 
the limits (7.1) without changing the variable x. As a result, we arrive at 

+ — ; = at, 



2(/xi - b)yja - JTi 2(/x 2 - b)yja - ~fT 2 (7 , o\ 



+ - — 2 = dx — b dt. 



2^Ja - JTi 2^/a - /j, 2 

'■•'"" ] 1 () " £7=Sj£ saving Dimpic puica ^ 2 , 



The latter equations include one differential of third kind Q2 = having simple poles Q 



Q2 and one differential of second kind Cl^ having a double pole at infinity on the Riemann surface 
P = {£ 2 = a — fj,}. Consider the mapping 

rPi rPi rPi rP2 

/ o«+/ nw=«i, / ft 2 +/ 2 = z 2 , ^ = (A*i,e*), (7-14) 

JPo JPo JPo JPo 

where, Z\ = /3 2 t + Z\q, z 2 = x + bt + z 2 o, 210, -220 =const and, as above, Po = («, 0). Integrating it 
explicitly, we obtain 

Inverting these relations yields the following formal solution for the HD equation 

U(x,t) = m + ii 2 = (6 + 6) 2 - 266 = *i - 2ft (ft - ^i) ^/ 2 ^_ 2 '/ 2 » (7.16) 
zi = x + 6t + zio, z 2 = fti + £20, zio, 2 20 = const. 

It is seen that U depends on x rationally (quadratically, as already mentioned above) and U is 
unlimited as x — ► ±00. 

However, this solution does not take into account the reflection phenomenon described above: 
when the variable /ii vanishes, the corresponding point Pi G IP jumps from the pole °f the 
differential f2 2 to Q^- According to mapping ( [7.14 ), this results in changing the phases z\,z 2 in 
( 7.16| ) by the constants 

01 ~ 02 



I 2 = 2ft, respectively q= [ * fi 2 = 2 log 



01+02' 

the latter being already defined in (|7.7| ). It follows that the actual solution U(x,t) to HD equation 
consists of an infinite number of pieces described by ( |7.16| ) with 

z x = x + bt + 2ft N + 210, z 2 = 2 t + Nq + z 20 , N G Z (7.17) 

and glued along peafc /ines {x = qisr(t)} in the plane (x,t), where for a fixed time t, the function 
QN(t) gives x-coordinate of iV-th peakon. Now if we assume Z20 to be imaginary and Ziq real, the 
series of pieces will provide a real bounded peakon solution. 

The functions (?jv(i) can be found as follows. Along the peak lines one of the variables n, say 
fix, vanishes implying dfi% = 0, 6 = 02- Putting this into ( 7.15| ), we find 

f^f = e*"« 1\ 02 + ^2 = Zl . 

6 + 02 
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Substituting here ( 7.17 ), putting x = qN(t), and eliminating £21 we obtain the sought expression 

1 + e Z2 - q / 2 



q N (t) = -bt - z 10 - 20iN + P1+P2 



1 



=22-9/2 ' 



Z2 = fat + Nq + Z20. 



It follows that for \t\ >> 1 any peak moves with constant velocity —b, and as t evolves from —00 
to 00 the peaks undergo the phase shift x — > x — 2/?2- 
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